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Abstract. The paper gives an accurate description of the qualitative structure of an admissible BV 
solution to a strictly hyperbolic, piecewise genuinely nonlinear system of conservation laws. We prove 
that there are a countable set which contains all interaction points and a family of countably many 
Lipschitz curves such that outside 3T U u is continuous, and along the curves in J?*, u has left 
and right limit except for points in 0. This extends the corresponding structural result in [8, 10] for 
admissible solutions. 

The proof is based on approximate wave-front tracking solutions and a proper selection of disconti- 
nuity curves in the approximate solutions, which converge to curves covering the discontinuities in the 
exact solution u. 

<N 

p 1 1. Introduction 

This paper is concerned with the qualitative structure of admissible solutions to the strictly hyperbolic 
^ N x N system of conservation laws in one space dimension 

+-> 

iu t + f(u) x = Hil+xR^SlcM" /eC 2 (fi,M), 
£j ' \u\t=o = u u Q eBY(R;n). 

CN| We assume the strict hyperbolicity in ft: the eigenvalues {Xi(u)}fL 1 of the Jacobin matrix A(u) — 

Df(u) satisfy 

Xi(u) < ■ ■ ■ < Ajv(u), well. 
Furthermore, as we only consider the solutions with small total variation, it is not restrictive to assume 
that f2 is bounded and there exist constants {Aj}jL , such that 

^ (1.2) A fc _! < A fc (u) < A fe , Vuefi, fc = l,...,iV. 

Let {ri(u)}fL 1 and {J^ (ii)}j!Li be a basis of right and left eigenvectors, depending smoothly on u, such 
that 

lj(u) ■ n(u) = Sij and \ri(u)\ = 1, i, j = 1, . . . , N. 
Let Ri[uo](uj) be the value at time oj of the solution to the Cauchy problem 

^=n(u(u)), u(0)=u o . 

We call the curve i?i[ M o] the i-rarefaction curve through uq. 

We say that the system (1.1) is piecewise genuinely nonlinear if the set where VAi • Vi — is covered 
by ki transversal manifolds: more precisely, 

Zi := {u € n : |VA 4 • r t {u)\ = 0} = Q z{, 

3=1 

where Zj is a N — l-dimensional manifolds such that 

(1) each Zf is transversal to the vector field r-j(w), i.e. 

(1.3) (V(VAj • r"j) • ri){u) 7^ foTueZ?; 
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(2) each rarefaction curve i?i[ito] crosses all the Zf, and moreover defining the points lo j [uq\ by 

i^M(u;> ]) e zj, 

then j i y lj j [uq] is strictly increasing. 
This implies that along R il A^ has a finite number of critical points. 
Denote with A^ the set of points u between Z\ and Zf +1 : 

A{ :={ue!]: uj 3 [u] < < oj j+1 [u}}. 
Without any loss of generality, we assume that 

(1.4a) VAi • ri{u) < if j is even, u 6 A^, 

(1.4b) VA* • n{u) > if j is odd, u e A{. 

From now on, we assume that every characteristic field of (1.1) is piecewise genuinely nonlinear. 

It is well known that, because of the nonlinear dependence of the characteristic speeds Aj(u) on the 
state variable it, the solution to (1.1) develops discontinuities within finite time, even with smooth initial 
data. Therefore, in order to construct solutions globally defined in time, one considers weak solutions 
interpreting the equation (1.1) in a distributional sense. We recall that u € C(M + ; L\ oc (&; R^)) is a weak 
solution to the Cauchy problem (1.1) if the initial condition is satisfied and, for any smooth function 
(j> € Cl(]0,T[xR) there holds 

cT r 

4>t{t, x)u(t, x) + 4> x (t, x)f(u(t, x))dxdt = 



u(t, x) 



As a consequence of the weak formulation, it follows that a function with a single jump discontinuity 

f u L if x < at, 
I « R if x > at, 

is a solution to (1.1), if and only if the left and right states it L ,it R G R , and the speed a satisfy the 
Rankine-Hugoniot condition 

(1.5) f(u K )-f(u h )=a(u R -u L ). 

By the strict hyperbolicity, it is known that for any u~ € f2 there exists so > and N smooth curves 
S^m - ] : [— S(b s o] — ^ associated with functions &i : [— Soi s o] ~~ * ^ such that 

(1.6) ai(8)[3i[u-](a) - u~] = /(fr^u"]^)) - f{vT) 
and satisfying 

Si[u-}(0) = u~ , &i(0) = Ai(u-), ^-5 i [«-](0)=r i («-). 

as 

The curve Si [it - ] the i-th Hugoniot curve issuing from u~~ and we also say that [it~,u + ] is an i- 
discontinuity with speed &i(u~ , u + ) := er(s) if ?i + = Si[tt — ](s). 

Since weak solutions to (1.1) may not be unique, an entropy criterion for admissibility is usually added 
to rule out nonphysical discontinuities. In [10], T.P. Liu proposed the following admissibility criterion 
valid for weak solutions to general systems of conservation laws 

We say the i-discontinuity [u~~ u + = Si[u~](s)], is Liu admissible if it satisfies Liu admissible 
condition: for sq > 

where u = Si[u~]{r) for each r €]0, Sq[, and for Sq < 

<Tj(u + , u~) > ai(u, u~), 

where u = Sj[u~](-r) for each r G]so, 0[. 

Let [u~,u + ], u + = Si[u~](s), be a Liu admissible i-discontinuity. Following the notation of [10], we 
call the jump [u~,u + ] simple if Vr G]0, s[ when s > (Vr g]s,0[ when s < 0), 

ai(Si[u~](T),u~) < ai(u,u~) (<Ti(S , l [u _ ](r), u~) > fri(u,u~)). 

If [u~ , u + \ is not simple, then we call it a composition of the waves [u~ , tti], [ui, u 2 ], ■ ■ ■ , [u;, it + ], if 

(1.7) u k = S l [u-](s k ) and d- t {u k , u^ 1 ) = a l {u + , u~), 
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where 



= Sq < Sl < S2 < • • • < si < s (or S < S; < • • • < Si < Sq = 0), 



and there are no other points r such that (1.7) holds. 

In [10], under assumption of piecewise genuinely nonlinearity, it is proved by using Glimm scheme 
that if the initial data has small total variation, there exists a weak BV solution of (1.1) satisfying Liu 
admissible condition. Therefore, it enjoys the usual regularity properties of BV function: u either is 
approximately continuous or has an approximate jump at each point (x, t) € M + x K \Af, where Af is a 
subset whose one-dimensional Hausdorff measure is zero. 

In [10], the author shows much stronger regularity that u holds. The set Af contains at most countably 
many points. Moreover, u is continuous (not just approximate continuous) outside Af and countably many 
Lipschitz continuous curves. 

In [8], the authors adopt wave- front tracking approximation to prove the similar result for (1.1) with 
the assumption that each characteristic field is genuinely nonlinear. Moreover, the authors where able to 
prove that outside the countable set Af there exist right and left limits u~, u + along the jump curves in 
the uniform norm, and these limits are stable w.r.t. wavefront approximate solutions: more precisely, for 
each jump point (not interaction point) of the solution, there exists a jump curve for the approximate 
solution converging to it and such that its left and right limit converge to u~ , and u + uniformly. In [6] 
(Theorem 10.4), the author generalize his result in [8] to the case when some characteristic field may be 
linearly degenerate. 

To prove this new regularity estimates one has to overtake additional difficulties, and this is the reason 
why they have so far been restricted to genuinely nonlinear of linearly degenerate systems: in fact the 
proof in [6] is based on the wave structure of the solution to genuinely nonlinear or linearly degenerate 
systems, where only one shock curve passes through the discontinuous point (which is not an interaction 
point) of the admissible solution. 

In this paper, we extend the techniques of [6] to prove an analogous result about global structure of 
admissible solution for piecewise genuinely nonlinear system of (1.1) by means of wave- front tracking 
approximation. This not only completes the corresponding result in [10], but also makes it possible to 
prove SBV regularity for the solution of piecewise genuinely nonlinear strictly hyperbolic system. In fact, 
one of the key argument for SBV regularity in the proofs contained in [4] and [5], is that outside the 
interaction points the left and right values of jumps are approximated uniformly by wavefront approximate 
solutions. 



As we said, the assumption of genuinely nonlinearity done in [8] implies that there is only one shock 
curve passing through the discontinuous point (which is not an interaction point) of the admissible 
solution. For piecewise genuinely nonlinear case, however, due the presence of the composite discontinuity, 
there may be several discontinuity curves passing through the same discontinuity point which is not an 
interaction point. For example, consider a scalar equation where / has two inflection points (it is thus 
clearly piecewise genuinely nonlinear), Figure 1, and let Uq be the initial data be 




Figure 1. 



Figure 2. 



u = 




if x < x\, 

if Xi < x < X2, 

if X2 < X < X3, 

if x > X3. 
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The same Figure 2 shows two shocks connecting value U\,u,2 and it3,it4 respectively interact with a 
center rarefaction wave and eventually have the same speed at P (which is thus not and interaction point) 
and combine together becoming a single shock. Clearly such wave pattern can not happen if / is convex 
or concave. 

In this paper we prove the following theorem: 

Theorem 1.1. Let u be a Liu admissible solution of the Cauchy problem (1.1). Then there exist a 
countable set O of interaction points and a countable family S? of Lipschitz continuous curves such that 
u is continuous outside and Graph(^). 

Moreover, suppose u(to, x) is discontinuous at x — Xq, and (to, xq) (fc 0. Write u L = u(to, Xq—), u r = 
u(to,xo+) and suppose that u R = Si[u L ](s) with s > (s < 0). 

• If [u L , u R ] is simple, there exists a Lipschitz curve y € 3 ' , s.t y(to) = xq 

u L = lim u(x,t), u R = lim u(x,t) 

Co;,*) — Ktch^o) (^i-H^o^o) 

and the curve y propagates with shock speed a(v}^,v^) at (to,xo), that is 

y(t )(u*-u h ) = f(u K )-f(u h ). 

and 

y(to) < &i(Si[u L ](T),u L ), Vr e M (y(t ) > ^(^[^(r), u L ), Vr e [s,0]). 

• If [ii L , u R ] is a composition of [u L , u%], [ui, U2], ■ ■ ■ ,[ui, u R ], then there exists p Lipschitz contin- 
uous curves j/i, • • ■ , y p € , p < I + 1 satisfying 

- yi (to) = ■■■ = y P (t ) = x Q , 

- y'i(to) = ■■■ = y' p (to), 

- yi(t) < • • • < y p (t),for all t in a neighborhood of to, 

s.t. 

u L = lim u(x,t), u R — lim u(x,t), 

x< B1 (t) x> Hp (t) 
(ic,t)->(t ,aiQ) (a,t)-»(t ,* ) 

and if in a small neighborhood of (to, xq), yj and yj+i are not identical, one has 
(1.8) Uj = lim u(x,t). 

(a=,t)->(*0' x 0) 

Also, these curves propagate with speed o^it 1 ",^) at (to,xo), that is 

y n (t )(u R - u L ) = f(u R ) - f(u L ), n € {1, • • • ,p}. 

and the stability conditions hold: 

y n (to) < 5- < (S i [u L ](r) 1 u L ) > Vr G [0,s], (y n (t Q ) > 5- i (#[u L ](T),u L ), Vr G [s,0]). 

As in [6], the above result is based on this strong convergence result for approximate wave-front 
solutions. 

Theorem 1.2. Consider a sequence of wave-front tracking approximate solutions u v (see Section 3 
for definitions) converging to u in L\ oc . Suppose P = (r, £) is a discontinuous point of u and write 
u L = m(t, £— ), u R = u(t, £+). Assume there are only I Lipschitz continuous curves 3~ 3 y n : [i^,f^"] 1— > 
K, n = 1, • • • , I passing through the point P and 

Ui(t) < • • • < yi(t) in a small neighborhood of V. 

Then up to a subsequence, there exists y n ^ : [i~ i—)- K, n = 1, • • • ,1 which are discontinuity curves 
of U v with uniformly large strengths, where t~ „ — > t~ , t+ v — > t+ and 

(t) -> y n (t) for every t € [t n , i+] . 

Moreover, one has 

lim limsup sup |wi,(£, t) — ti L | ] =0, 

r ~ ! '°+ !/— >oo I as<»l,i/C*) 
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lim limsup sup \u v (x,t) — u \ =0. 

r^-0+ ,,-j.oo I „>y, il ,(i) J 

A brief outline of this paper follows. 

In Section 2, we recall the definition of construction of Riemann solvers introduced in [3]. 

In Section 3, we briefly describe the wave-front tracking approximate scheme which originally designed 
for general strictly hyperbolic system (sec [1]). In particular, we introduce the definition of interaction 
and cancelation measures. 

Section 4 contains the main idea of the paper: the definition of subdiscontinuity curves and (e, k)- 
approximate subdiscontinuity curves in the approximate wavefront solution. In this section we show that 
their number is uniformly bounded with respect to approximation parameter. 

In Section 5, we finally give the proof for Theorem 1.1 and Theorem 1.2. 

In Section 6, we construct a strictly hyperbolic 2x2 system of conservation laws, which is not piecewise 
genuinely nonlinear and its admissible solutions to some initial datum do not have the structural properties 
described in Theorem 1.1. 



2. Solution of Riemann problem 

As in [2], for a fixed point u° E il and i £ {1, ■ ■ ■ ,N}, one can construct smooth vector- value maps 
fi = fj(u, Vi, Oi) for [u, Vi, <Ti) E M. N xKx! with fi(u, 0, a) = ri(u) for all u, Ui. Setting l\ := h(u°), we 
can normalize fi such that 



(2.1) «i,o-i)) 



1 i = j, 
i^j. 



Writing the speed function := If ■ D f (u)fi(u,Vi, <Ji) , we consider the set, for some So, Co > fixed 
and s > 

T l (s,u-) :={ 7 e LipaO,*],^ 2 ), 7 (0 = (u({)MZ)Mt)) 
u(0)=u-,|tt(T)-«-|=r,« i (0) = Q J 
\v i \<S 1 ,\a i {r)-X i (u )\<2C S i y 
Given a curve 7 € 1^, we define the scalar flux function 

(2-2) Mr, 7) = r \i(u(£)M£)M0)<%- 



Moreover, we define the lower convex envelope of fi on [a, b] C [0, s] as 
cony fi(r; 7) := inf {#/,(/, 7) + (1 + 9) for", 7); 

[a,b\ l 

9e [0,1], t',t" € [a,b],T = 0T / + {1 + 6)t"\. 

Then define a nonlinear operator % tS : Ti(s,u~) — > Ti(s,u~) by setting 7 = 7:= (fa, V{, (x,),where 

u(t) =u~ + J q t fi(u(0, Vi(Q, cr(0)d^ 
(2-3) < Ui(r) = /i(r; 7) - conv [0iS ] f l (r, 7) , 

= ^ conv [o, s ] /i( r ;7)- 

One can show that 7" is a contraction in Ti(s,u~) with respect to the distance 
where 

7 = (u,^,^), 7' = (u',-y-,cr-) e r,(s,?0- 
Hence, for any s and u~ in a small neighborhood of u , % has a unique fixed point, which is a Lipschitz 
continuous curve 

7(t) = (u(T;u~,s),Vi(T;u~,s),ai(T;u~,s)), r e [0,s]. 
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Then the elementary curve for z-th family is defined as 
(2.4) Ti[u-](s):=u(s;u-,s). 
After adopting the following notations 

CTiki~](s,T) : = ai(r;u~ ,s), 

(2-5) /i[«-](«,T):=/i(r;7). 

and recalling that the Riemann problem is the Cauchy problem (1.1) with piecewise constant initial data 
of the form 

(2-6) u ( a! ) = | mR x>Q; 

where m l , zz R are two constants. One has the following theorem [3]. 

Theorem 2.1. For every u £ Q and s > sufficiently small, 

(1) N Lipschitz continuous curves s >— ¥ Ti[u](s) £ SI, i = 1,...,N, satisfying lim s _>o = 

(2) jV Lipschitz continuous functions (s, r) h- > <7,[u](s, r), with < r < s and i = 1, . . . , N , satisfying 
t i— > <Tj[u](s, t) which are increasing and such that <Ji[u](s, 0) = Aj(ii), 

wz'zTi i/ie following properties. 

When u L £ fl, u R = Ti[u L ](s), for some s sufficiently small, the unique Liu admissible solution of the 
Riemann problem (l.l)-(2.6) is defined a.e. by 

!u L x/t < CTi[u L ](s, 0), 

Ti [u L ] (r) x/t = o-i [u L ] (s, t) , r £ [0, s] , 
u R x/i > (Ti[u L ](s, s). 



For the case when s < 0, a right state w R = T;[u L ](s) can be constructed in the same way as before, 
except that one replaces convr 0iS ] fi in (2.3) with the upper concave envelope of /j on [s, 0]: 

cone 7,(1-5 7) := sup f Of if/, 7) + (l+8)/ ( (/,7); 

[0,1],t',t"g [a,6],r = ^' + (1 + ^)/'}, 

and looks at the fixed point of of the integral system (2.3) on the interval [s, 0]. 

Because of the assumption (2.1) and the definition (2.4), the elementary curve T, [u L ] is parameterized 
by its z-th component relative to the basis Tx(u ), ■ ■ ■ ,rj^(u ) i.e. 

(2.8) 8 = <I»,T < [u L ](*)-« L ). 

Remark 2.2. In [3], it is proved that if u L ,ii R £ A^" with some k even (odd) and u R = Ti[u L ](s), s > 
(s < 0), the solution u of the Riemann problem with the initial date (2.6) is a center rarefaction wave, 
that is for t > 0, 

!u L if x/t < X x (u L ), 

u R tfx/t>\ l (u R ), 
Ri[u h }(r) ii x/t £ [Xi(u L ),Xi(u R )], x/t = X i (R i [u L ](T)), 

where r £ [0, s] (r £ [s, 0]) such that s = (Z°, i?i[u L ](s) — u L ). Notice that u is smooth for t > 0. 

Remark 2.3. As shown in [3] (also see Remark 4 in [1] and Section 4 of [10]), under the assumption of 
piecewise genuine nonlinearity, the solution of the Riemann problem provided by (2.7) is a composed 
wave of the z-th family containing a finite number of rarefaction waves and Liu admissible discontinuities. 
Recalling Theorem 2.1, one knows that the regions where the fi-component of the solution to (2.3) 
vanishes correspond to rarefaction waves, while the regions where the Wi-component of the solution to 
(2.3) is different from zero correspond to admissible discontinuities. 
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The Liu admissible solution [3] of a Riemann problem for (l.l)-(2.6) is obtained by constructing a 
Lipschitz continuous map 

(2.9) s := ( Sl , . . . , s N ) h> T[u L ](s) := T N [T N . X [■■■ [Ti[u L ](«i)] • • • ] {a N -i)] = « R , 

which is one to one from a neighborhood of the origin onto a neighborhood of u L . Then we can uniquely 
determine intermediate states u L = u)q, u>i, . . . , u>n = u R , and the wave strength s\, S2, ■ ■ ■ , sjv such 
that 

Wi=Ti[u)i-i](si), i = l,...,N, 

provided that \u L — u K \ is sufficiently small. 

By Theorem 2.1, each Riemann problem with initial date 



(2.10) u 



Ui-i x < 0, 
uji x > 0, 



admits a self-similar solution u,-, containing only i- waves. We call the i-th elementary composite wave 
or simply i-wave. Therefore, under the strict hyperbolicity assumption, the solution of the Riemann 
problem with the initial data (2.6) is obtained by piecing together the self-similar solutions of the Riemann 
problems given by (1.1)-(2.10). 

Indeed, from the strict hyperbolicity assumption (1.2), the speed of each elementary i-th wave in the 
solution Ui is inside the interval [Aj_i, Aj] if s <C 1, so that the solution of the general Riemann problem 
(l.l)-(2.6) is then given by 

u L x/t < A , 

Ui(x,t) Xi-x < x/t < Xi,i= 1,...,N, 
u R x/t > Aat. 



(2.11) u(x,t) 



3. Description of wave-front tracking approximation 

In [1], the authors provide an algorithm of wave- front tracking approximation for vanishing viscosity 
BV solutions to the strictly hyperbolic system which is much more general than the case discussed here. 
We modify a little about its algorithm in order to simplify our analysis. Due to Theorem 2.1, one knows 
that the solution constructed by such approximation is Liu admissible. 

Wave-front tracking approximation is an algorithm which produces piecewise constant approximate 
solutions to the Cauchy problem (1.1). In order to construct approximate wave-front tracking solutions, 
given a fixed e > 0, we first choose a piecewise constant function Uq which is a good approximation to 
initial data uq such that 

(3.1) Tot.Var.{w } < Tot.Var.{u }, \\u e - u \\ L i < e, 

and Uq only has finite jumps. Let x\ < • • • < x m be the jump points of Ug. For each a = 1, . . . , m, we 
approximately solve the Riemann problem (just shifting the center from (0, 0) to (0, x a )) with the initial 
data of the jump [u (x a —), [u : (a; Q ,+)] by a function w(x,t) = </>( ) where is a piecewise constant 
function. The straight lines where the discontinuities locate are called wave- fronts (or just fronts for 
short). The wave-fronts can prolong until they interact with other fronts, then at the interaction point, the 
corresponding Riemann problem is approximately solved and several new fronts are generated forward. 
Then one tracks the wave-fronts until they interact with other wave-fronts, etc... In order to avoid the 
algorithm to produce infinite many wave-fronts in finite time, different kinds of approximate Riemann 
solvers should be introduced. 

3.0.1. The approximate i-th elementary wave. Suppose that itj is an i-th elementary composite wave 
which is obtain by solving Riemann problem with initial data (2.10) where u>i = Tj[tai_i](sj). For 
notational convenience, we write crj(r) := ai[uii-i]{si, r). Let 

o-i(si) - o~i (0) 



P 



1 



and 



o-i(0) + - [o-i(si) - o-i(Q)\ , I = 0. • • • ,p - 1. 
P 
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We set 

Ui-lJ = Ti[u)i-i](8i t t), 

where 

_ \mm{s € [0,Si],ai(s) =$i}, > 0, 
| max{s G [sj,0],(Tj(s) = < 0. 

Then the i-th elementary composite wave Ui is approximated by Ui as the following, 

{Ui-i x/t < i? ijQ , 

Ui-x,u ■d i -x,i-x<x/t<ti it i, (I = !,■■■ ,P-1), 
cjj ar/* > i?i. p -i. 

Notice that Ui consists of p fronts with small strength. 

3.0.2. Approximate Riemann solver. Suppose at the point (tx,xx), a wave-front [u L ,u M ] of strength s' 
belonging to i'-th family interacts from the left with a wave-front [u M ,u R ] of strength s" belonging to 
i"-th family for some i', i" G {1, • • • , N} such that 

u Af =T i /[u L ](« / ), v R = T i „[u M ]{s"). 

Assume that \u L — u R \ sufficiently small. Then at the interaction point, the Riemann problem with the 
initial data of the jump [ii L ,u R ] may be solved by two kinds approximate Riemann solver according to 
different situation. 

• Accurate Riemann solver: It replaces each elementary composite wave of the exact Riemann 
solution (refers to Ui in the solution (2.11)) with an approximate i-th elementary wave defined 
by (3.2). 

• Simplified Riemann solver: It only generates the approximate elementary waves belong to i'-th. 
and i"-th families with the corresponding strength s' and s" as the incoming ones if i' ^ i" or 
the approximate i'-th. elementary waves of strength s' + s" if i' = i". The simplified Riemann 
solver collects the remaining new waves into a single nonphysical front, traveling with a constant 
speed A, strictly larger than all characteristic speeds. Therefore, usually the simplified Riemann 
solver generates less outgoing fronts after an interaction than the accurate Riemann solver. 

Since the simplified Riemann solver produces nonphysical wave-fronts and they can not interact with 
each other, one needs an approximate Riemann solver defined for the interaction between, for example, 
a physical front of the i-th family with strength s, connecting u M , u R and a nonphysical front (coming 
from the left) connecting the left value u h and u M traveling with speed A. 

• Crude Riemann solver: It generates an approximate i-th elementary wave connecting u L and 
u M = Ti[u L ](s) and a nonphysical wave-front joining u M and u R , traveling with speed A. In 
the following, for simplicity, we just say that the non-physical fronts belong to the (N + l)-th 
characteristic field. 

Remark 3.1. It is not restrictive to assume that at each time t > 0, at most one interaction takes place, 
involving exactly two incoming fronts, because one can always slightly change the speeds of the incoming 
fronts if more than two fronts meet at the same point. It is sufficient to require that the error vanishes 
when the approximation solutions converge to the exact solution. Actually, suppose x=y(t) is a front in 
an approximate solution u with parameter e and u L — u(t, x—) and u R = u(t, x+) such that 

(3.3) u R = T i [u L ](s) 

for some index i G {1, • • • , N} and wave strength s. Then the following holds 

(3.4) \y(t)-o-i[u L ](s,T)\ <2e, VtG[0,s]. 

Remark 3.2. There are three kinds of physical wave-fronts. Suppose [it L ,M R ] is wave front in an approx- 
imate solution, and it R = Tj[u L ](s) (s > 0). Recalling the notaton (2.5) and Remark 2.3, and writing 
fi(r) :— /i[it L ](s, t), one knows that [u L ,u R ] may be one of the following three kinds of fronts: 

• Discontinuity front if fi(r) > conv]o, s [ /i( r )i ^ r G]0i s[\S, where S is set of finite cardinality. 
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• Rarefaction front if fi{r) = convio jS r fi(r), Vt €]0,s[. In this case, u h ,u R G Af for some k even 
and m r = Ri[u L ](s). 

• Mixed front if there exist ]a, b[C]0, s[ such that Ji(r) > convio, s [ fifr), Vr e]a, b[ and ]c, d[C]0, s[ 
such that /j(r) = conv ]0:S[ / ( (r), Vr e]c,d[. 

3.0.3. Interaction potential and BV estimates. In order to check the total variations of approximate 
solutions are uniformly bounded with respect to time, one needs the estimate on the difference between 
the strength of the incoming waves and the strength of the outgoing waves produced by an interaction. 
Suppose two wave- fronts with strength s' and s" interact and fl, f" are the corresponding scalar flux 
function defined by (2.2). We define the amount of interaction X(s' ,s") between s' and s". 
When s' and s" belong to different characteristic families, set 

(3.5) l(s',s") = \s's"\. 

When s' , s" belong to the same family, 

(a) If s" > 0, we set 

l(s',s") := f | covvfXt) - r conv (fi U 

S ^ | conv(/V)/f (£ - s')) - conv {][ U f{ )(0\d£, 

i [0,s"J [0,s'+s"\ 

(b) if -s' < s" < 0, we set 
Hs',s") := f 5 | conv £(£) - conv 

JO [0,S>] [0,8' +8"] 

I 

+ [ S | cony cone 

Js'+s" [°> s 1 [a'+a",a'] 

(c) if ,s" < — s', we set 
I( S ', S ") := f S | conc/r(0 - cone 

Js" W'fi] [a" -a'] 

|conv(/'(0)- conv / w (0|C 

[a",0] [-a',0] 

Throughout the paper, we write A < B (A > B) if there exists a constant C > which only depends 
on the system (1.1) such that A < CB (A > CB). 

Recall the Lipschitz continuous map T defined in (2.9) and suppose u M = Ti[u L ](si), u R = Tj[u M ](s2) 
and u R = T[u L ](s). By Glimm's interaction estimates proved in [2] (also see Lemma 1 in [1]), one has 

(3.6) |s-si -s 2 | <X(si,s 2 ) 

where = (si, S2, • • • , %), h = 1,2 with s n — for n =/= h. 

At each time t > when no interaction occurs, and the approximate solution u has jumps at X\, . . . , x m , 
we denote by 

U)\ , . . . , LO m , S\ , . . . , S m , %\ , . . . , if/I, 

their left states, signed strengths and characteristic families respectively: the sign of s a is given by the 
respective orientation of dTi[u]{s)/ds and r,;, if the jump at x a belongs to the i-th family. The total 
variation of u v will be computed as 

ot 

Since Tj[u°] is a Lipschitz continuous function and the Lipschitz constant is uniformly bounded for 
any ti° gO, one has 

(3.7) Tot.V&v{u(;t)}<V(t). 
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Then the estimate of increasing of Tot.Var{it„(-, £)} turns out to be the estimate of the total amount of 
interaction. Following [2], we define the Glimm wave interaction potential as follows: 

2( f ) : = Z_, I SqS ' 3 I + 4 z2 \a if! [ujp](sp,T") - a ia [uj a ](s a ,T')\dT'dT" . 

*a>*f> i Q=4f) <Ar+l"' JO 

Denoting the time jumps of the total variation and the Glimm potential as 

AV(r) = V(t+) - V(t-), AQ(r) = Q(r+) - Q(r-), 

the fundamental estimates are the following (Lemma 5 in [1]): in fact, when two wave- fronts with strength 
s' , s" interact, 

(3.8a) AQ(t) <Z(s',s"), 



(3.8b) A7(r) <l(s',s"). 

Thus one defines the Glimm functional 

T(t) := V(t) + CoQ(t) 

with Co suitable constant, so that T decreases at any interaction. Using this functional, one can prove 
that their total variations are uniformly bounded. Moreover, one can also show that the number of 
wave-fronts remains finite for all time (see section 6.1 of [1]). This makes sense for the construction of 
approximate wave-front tracking solutions. 



3.0.4. Construction of the approximate solutions and their convergence to exact solution. The construc- 
tion starts at initial time t = with a given e > 0, by taking ito !£ as a suitable piecewise constant 
approximation of initial data ito, satisfying (3.1). At the jump points of Mo, £ , we locally solve the Rie- 
mann problem by accurate Ricmann solver. The approximate solution then can be prolonged until a 
first time t\ when two wave-fronts interact. Again we solve the Riemann problem at the interaction 
point by an approximate Riemann solver. Whenever the amount of interaction (see Section 3.0.3 for the 
definition) of the incoming waves is larger than some threshold parameter p = p(e) > 0, we shall adopt 
the accurate Riemann solver. Instead, in the case where the amount of interaction of the incoming waves 
is less than p, we shall adopt the simplified Riemann solvers. The threshold p is suitably chosen so that 
the number of wave-fronts remains finite for all times. And we will apply the crude Riemann solver if 
one of the incoming wave- front is non-physical front. One can show that the number of wave fronts is 
uniformly bounded (see Section 6.2 in [1]). 

We call such approximate solutions e-approximate front tracking solutions. At each time t when there 
is no interaction, the restriction u e (t) is a step function whose jumps are located along straight lines in 
the (x, i)-planc. 

Let {ci/}^! be a sequence of positive real numbers converging to zero. Consider a corresponding 
sequence of e„-approximate front tracking solutions u v :— u tv of (1.1): it is standard to show that 
the functions t M- u„(t, •) are uniformly Lipschitz continuous in L 1 norm, since (3.7) and (3.8) hold 
independent of the parameter e u , u u (t, •) have uniformly bounded total variation. Therefore by Helly's 
theorem, u„ converges up to a subsequence in L 1 1 oc (i? + x IR) to some function u, which is a weak solution 
of (1.1). 

It can be shown that by the choice of the Riemann solver in Theorem 2.1, the solution obtained by 
the front tracking approximation coincides with the unique vanishing viscosity solution [3]. Furthermore, 
there exists a closed domain D C i 1 (K, ft) and a unique distributional solution u, which is a Lipschitz 
semigroup T> x [0, +oo[— )> T> and which for piecewise constant initial data coincides, for a small time, 
with the solution of the Cauchy problem obtained piecing together the standard entropy solutions of the 
Riemann problems. Moreover, it lives in the space of BV functions. 

For simplicity, the pointwise value of u is its L 1 representative such that the restriction map t > u(t) 
is continuous form the right in L 1 and x i-> u(x,t) is right continuous from the right. 
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3.0.5. Further estimates. To each u u , we define the measure ij} v of interaction and the measure ^t* of 
interaction and cancelation concentrated on the set of interaction points as follows. If two physical fronts 
belonging to the families € {1, . . . ,N} with strength s' , s" interact at point P, we denote 

(3.9a) nl({P}):=l(s',s"), 

(3.9b) tf?({P}) ■■=!(*',*")+ | I s ' 1 .„ 

The wave strength estimates (3.6) yields balance principles for the wave strength of approximate 
solutions. More precisely, given a polygonal region T with edges transversal to the waves it encounters. 
Denote by W* in , W^ ut the positive (+) or negative (— ) i-waves in u v entering or exiting T, and let 
Wl in = Wl\ n — W^~ ln , W* out = W*+ Ut — Wl~ out . Then the measure of interaction and the measure of 
interaction-cancelation control the difference between the amount of exiting i-waves and the amount of 
entering «-waves w.r.t. the region as follows: 

l^ut-^l<Ml c (r). 

The above estimates are fairly easy consequences of the interaction estimates (3.8) and the definition of 
A 4 !/) Av • On the other hand, the uniform boundedness of Tot.Var.{w(-, t)} w.r.t. time t and parameter v 
implies that ^ and fi — v are bounded measures for all v 

By taking a subsequence and using the weak compactness of bounded measures, there exist bounded 
measures fj, 1 and /i IC on E + x E such that the following weak convergence holds: 

4. Construction of subdiscontinuity curves 

Suppose [u L , u R ] with u R = Ti[u h ]{s) is a wave front of i-th family in the approximate solution u v , and 
the wave curve r i-> Ti[u L ](T) (see Theorem 2.1) intersects Zj,--- , Zf +P at Uj, ■ ■ ■ ,Uj +p for < r < s, 
such that 

Uj = T l [u L ](s i ). 

Then we say that the wave front [u L ,u R ] has (i, fc)-substrength sf := Sk+i — Sk and we decompose the 
front into (i, fc)-subdiscontinuity fronts with strength s^, where fee {j, ■ ■ ■ ,j+p—l}<~)2Z when s > 0, 
or k € { j, ■ ■ ■ , j + p — 1} n (2Z + 1) when s < 0. The points and Ufc+i are connected by the part of 
the curve Ti[u L ](-) inside Af . 

We denote the family of all (i, /c)-subdiscontinuity fronts as . 

It is obviously that only mixed fronts and discontinuity fronts can have (i, A:)-substrength > for 
some k, which means they can be decomposed into subdiscontinuity fronts. 

Lemma 4.1. In a wave-front tracking approximation solution, an interaction can only generate at most 
one subdiscontinuity front with strength s^ for some i, k. 

Proof. By the construction of approximate Riemann solver and the uniformly small total variation of 
approximate solutions, it is sufficient to prove that the Lipschitz continuous curve Tj[u°](-) : [0, s] —> R N 
can intersect with Zf at most once for any «° £ SI and all j if s > is sufficiently small. In fact, by 
Theorem 2.1 lim t ^ Ti[u°](t) = r^u ), one has for t e [0, s], 



St. 



Recall the assumption Vi(u) ■ n > on Zf, one has 

jTi [u°] (t) ■ n > on Z{ for t G [0, s] 
as long as s small enough. This concludes that Tj[u°](-) : [0, s] —> R N can intersect Z\ at most once. 

□ 
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Suppose 2/1 G , y2 G <Sf and j/i, J/2 belong to the same wave front, then we say artificially y\ is on 
the left (right) of y 2 if k',k" are even (odd) and k' < k" (k' > k"). Then it is easy to see the following 
lemma. 

Lemma 4.2. Suppose two subdis continuity fronts y\ G <S* , yi G S\ interact and generate two subdis- 
continuity fronts y[ G S!f , y' 2 G 5* . TTien i/fc' < fe" (fc' > k"), y\ must be on the left (right) of y 2 and 
also y[ must be on the left (right) of y' 2 . 

Base on these two lemmas, we can follow the idea of [7] to define approximate subdiscontinuity curves. 

Definition 4.3. Given e ^ 0, in u v , an (e, i, fc)-approximate subdiscontinuity curve is a polygonal line 
in (x, t)-plane with nodes (to,Xo), (ti,xi), • • • , (t n ,x n ) satisfying 

(1) (tj, Xj) are interaction points with < to < t\ < ■ ■ ■ < t n . 

(2) For 1 < j < n the segment joining (tj—i,Xj—i), (tj,Xj) is an (i, fc)-subdiscontinuity front with 

> e/2 when e > (sf < e/2 when e < 0), and there is at least one index j' G {1, ■ • • , n} 
such that the wave front connecting (tj'-±,Xf-i), (tj',Xj/) has (i, A:)-substrength > e if e > 
(sf < e if e < 0). 

(3) Vfc < AT, one selects the (i, /c)-subdiscontinuity fronts with larger speed, that is, if two fronts with 
substrengths > e/2 interact at the node (xk,tk), then the front of the (e, i, fc)-approximate 
subdiscontinuity curve is the one coming from the left. 

An (e, z, fc)-approximate subdiscontinuity curve which is maximal w.r.t. set inclusion is called a maximal 
(e, i, k)- approximate subdiscontinuity curve. 

Let Mf v (e) be the number of maximal (e, i, fc)-approximate subdiscontinuity curves in u v . 

Lemma 4.4. For fixed k and e, Mf v {e) is uniformly bounded w.r.t v. 

Proof. Wc consider the case for e > (and the case when e < is similar). Since the total variation 
of u v (0, •) is uniformly bounded by (3.1), the number of maximal (e, i, fc)-approximatc subdiscontinuity 
curves which start at time t = is clearly of order e . The number of (e, i, /^-approximate subdisconti- 
nuity curves which start at a time to > and do not end in finite time is also of order e _1 , because the 
total variation of u v (,t) is uniformly bounded w.r.t time t and v. 

Now considering an (e, i, fc)-approximate subdiscontinuity curve y v which starts at a time to > and 
ends in finite time, we claim that 

(4-1) tf( Vv ) > 6 2 . 

If the claim is true, as the total amount of interaction and cancelation in the solution u v is uniformly 
bounded, the number of such (e, i, fc)-approximate subdiscontinuity curves is of order e~ 2 . 
Thus, combining these situations, we finally obtain the estimate 

which is uniformly valid as v — > oo. 

Now we prove the claim. Suppose one (i, fc)-subdiscontinuity front on y v with {i, fc)-substrength a > 
interact with a front of j'-th family with strength /3* at point P, generating an i-wave with (i, k)- 
substrength 7 > which means that either there is a front of z-wave with (i, &)-substrength 7 > or 
there is no fronts of z'-wave having (i, /c)-substrength. We also assume that 7 < a. Setting 9 = a — 7. 

First, we consider the case when i ^ j, we assume that i > j, the case i < j is similar to prove. For 
notational convenience, we also denote a, (3* , 7 as the front themselves. 

Assume that a locates on the front [u L ,u M ] with u M — T^u^^a*) and u R = Tj[u M ](/3*), 7 locates on 
the front of the i-wave [u L ,u R ] with u R = Tj[w L ](7*). 

We know that (see section 9.9 of [9]) 

(4.2a) u L = u L + 8j r M M ) + 0{\s L \)a* + o(\s L \), 

j<i 

(4.2b) u R = u R - Sj rj(u M ) + 0(\s R \)(3* +o{\s R \), 

3>i 
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where 

s L = (si,--- ,Si_i,0, ••• ,0) and s R : = (0, • • • , 0, s i+1 , ■ ■ ■ , s N ). 
Then, by (3.6), the assumption i > j and \u M — u R \ < |/3*|, we have 

(4.3a) \u L -u L \<\p*\+T{a\n, 

(4.3b) |u M -ii fl | < + Z(a*,/3*). 

From Glimm's interaction estimates (3.6), the parametrization (2.8) and the estimates (4.3), one con- 
cludes that the difference of (i, fc)-substrength between a and 7 is controlled by the amount of interaction 
and the strength of the wave /3* , that is 

(4.4) 6 = a- 7 <X(a*,P*) + \p*\. 

From definition of X (see Section 3.0.3), we know that here X(a*,/3*) = |a*/3*|. So we get |/3*| > 9 
since \a*\ <C 1. Therefore from |cc*| > e/2, by the estimates in (i)-(iii), we obtain 

T(a*,0*) > e0. 

By notation of interaction measure (3.9a), one obtains 

(4-5) rf,{P) > e9. 

Next we consider the case when i = j and a is on the left of /3*. (It is similar for the case when a is 
on the right of /?*.) 

First we assume that a > 0, /?* < 0. Since the decreasing of the (i, A;)-substrength can only be caused 
by interaction and cancellation effects, then similarly one has 

Second, we assume that a > 0, j3* > 0. 
From the equality (4.2a), one know that 

(4.6) (l a ,ii L -u L ) <I(a*,/3*). 
From the equality (4.2b), 

(l°,u M -u R ) + {l°,u M -u R ) 
= (l°,u R -u R ) <I(a*,/3*). 

Since (l°,u M - u R ) = -p* < 0, one has 

(4.7) (l°,u M -u R ) <2(a*, /3*). 

Noticing that (4.6) and (4.7) implies that the decreasing of the (i, fc)-substrength can be controlled by 
the amount of interaction, one has 

6<X(a*,/3*). 

Therefore from the definition of (3.9b), in the case when i = j one has 
(4-8) ^ C {P) > 9. 

From (2) of Definition (4.3), the (i, fc)-substrength all front outgoing from the terminal point of y v 
must be less than e/2 and there is at least one front on y v has (i, fc)-substrength larger than e. Then by 

(4.5) and (4.8), one can conclude that the claim (4.1) is true. □ 

Up to a subsequence, one can assume that M^ u {e) = M^(e) is a constant independent of v . Denote 

V k m % ■■ [*££7 , t^+] -> R, m = 1, • • • , M* (e) 

as the maximal (e, i, A;)-approximate subdiscontinuity curves in u v . 

Define =^„(e) as the collection of all maximal (e, i, /c)-approximate subdiscontinuity curves in u v for 
fixed e, i and k, i.e. 

^) = {vtu- m = !,■■■ ,Mf(e)}. 
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Set 

fc,e 

as the family of all approximate subdiscontinuity curves of i-th family in u v . 

Up to a diagonal argument and by a suitable labeling of the curves, one can assume that for each fixed 
h, m, as v — > oo the Lipschitz continuous curves y^j e „ converge uniformly to some Lipschitz continuous 
curves y^f which are called (e,i, k)- subdiscontinuity curves. Let us denote by 

^{e):={y k ^: m = 1, • • • M? (e)} 

the collection of all these limiting curves for fixed i, k, e. 
Let 

ft:=]J^ k (e). 

k,e 

denote the collection of all these i- subdiscontinuity curves. 

Lemma 4.5. Let y^ :]t~,t+[— > K be an (e,i, k) -subdiscontinuity curve. If t is such that 

(t, ^ ®j i/ien the derivative y„(t) exists. 

Proof. There exist eo > and y^ n v £ ^^(eo), such that 

as v — y oo and the substrength of y 1 ^ v is \s k \ > £q. 

Since in the wave-front tracking approximation, the change of speed of z-subdiscontinuity fronts is 
controlled by the measure ^1° ' . 

Then for any <5„ — > 0, we deduce that 

limsup sup \vi, tU (t) - y^ v {t')\ = 0. 

v-yoo |t— i'|<4„ 

From uniformly convergence y 1 ^ v — > y^ on a neighborhood of t 1 we obtain 
(4-9) y k m (t) = Km y k m , v . 

□ 

Recall the definition of generalized characteristics which will be used in the proof of Theorem 1.1. 

Definition 4.6. A generalized i- characteristic associated with the approximate solution u v , on the time 
interval [ti,^] C [0, oo), is a Lipschitz continuous function \ '■ [^1,^2] ~~ * (—00,00) which satisfies the 
differential conclusion 

X(t) € [\i(u„(x+,t)), X^Uvix-^))}. 

For any given (T, x) £ R, we consider the minimal (maximal) generalized i- characteristic through 
(T, x) defined as 

X~ {t) = min(max){x(t) : x is a generalized characteristic, x{T) = x}. 

The properties of approximate solutions yield that there is no wave-front of z-th family crossing x + 
from the left or crossing x~ from the right. 

Suppose ,%, u 3 y' v : [f~,t'+] -> K and ,%, v 3 y'l(-) : [t"~,t'l+] -> R. By Lemma 4.2 and the definition 
of (e,z, fc)-approximate subdiscontinuity curves, it turns out that either 

y' v (t)< v'M Vte[i*-,i>+]n[%-X + l 

or 

y' v {t)> y'l(t), Vte[t>-,t'+]n[t'r,€ + }- 

This makes the following definition well defined. 

Definition 4.7. Suppose y' £ 2?^ (e'), y" £ ST k (e") where k' and k" are both odd numbers (or even 
numbers), and y' v — > y', y" — > y" as v — > 00 and assume there exists a point (ioi^o) with to > such 
that 2/ (to) = y"(to) — xq, we say y' -< y" if there exists a neighborhood {t~ , t + ] of time to such that 
• y'(t) < y"{t) for all t £ [t~,t + ], 
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• either there exist t* G [t ,t + ] such that y'(t*) < y"(t*) 

or for all t G [t~,t + ], y'(t) = y'(t), and ki < k 2 (k\ > k 2 ). 

The next lemma rules out the situation when two subdiscontinuity with different sign of strength is 
tangent at the same point which is not the atom point of interaction and cancellation measure. 

Lemma 4.8. Suppose that y\ G i^ fel , J/2 G S?^ 2 with k\ is even and k 2 is odd, and assume that y± and 
7/2 pass through the same point (xo,to) with (j, ic ({(xq, to)}) — and there is no subdiscontinuity curve j/o 
such that 

Vi(t) < y {t) < yi (t) 
for any neigbourhood of to. Then j/i(io) ^ J/2 (£())• 

Proof. Suppose i^^ 1 3 y\. v — ¥ y\, Sf^l 3 y 2iL , J/2 and ti <v ,t2,v — > to- Let us denote the points 
yi,v(ti tV ), y 2 .vit 2M ) as A V ,B V respectively. Since there is no subdiscontinuity curves between y\ and y 2 , 
the strengths of all fronts crossing the segment A V B V tend to zero. 

Moreover, the total strength of the fronts of the other family tends to zero. In fact, if not, either they 
are canceled in the neighborhood of (x ,t ) or interacted with y\, v {t\, v ) ■> y 2l v(t 2 ^), which implies the 
uniform positivity of [i„ on a small region r„. This contradicts the assumption that /j, ic ({(x^, t )}) = 0. 
Therefore, the values of each u v along the segment A V B V remain arbitrary close to z-rarefaction curve. 

On the other hand, since the sign of the strengths of j/j and y 2 are different, one can always 
find A' vl B' v on A V B V and a positive constant c such that u u (A' u ), u v (B' v ) G for some k and \u v {A' v ) — 
u v {B' v )\ > c, therefore 

\\i(u v {A v )) - Ai(it„(^))| > c. 
Up to subsequence, we can assume that for all v 
(4.10a) \i{u v {A' v )) - \{u v {B' v )) > c/2, 

(4.10b) or Xi(uy(A' v )) - \i(u„(Bl)) < c/2. 

Let us consider the case (4.10a), the other case is analogous to prove. We take \ + through A' u , \ ~ 
through B' v . Since A'^ and B' v are in the same A^, if no uniformly large interaction occur on x + , \~ , 
they will interact with each other. We consider the region r„ bounded by A v B Vl \ + and \ ~ . Since no 
fronts can leave Y v through x + and X - By (3-5) and (3.9), we obtain that /i 7 (r i/ ) > c which contradicts 
the assumption fJ. IC ({(x , t )}) =0. □ 

5. Proof of Theorem 1.1 

Before proving the theorem, we recall the definition space-like curve. 

Definition 5.1. Let A be a constant larger than the absolute value of all characteristic speed. We say a 
curve x — y(t), t G [a, b] is space-like if 

\y{t 2 ) - y{ti)\ > \{t 2 - tx) for all a < t x < t 2 < b. 

From the definition one knows that any fronts can cross a space-like curve at most once. 

Proof of Theorem 1.1. Let O consists of all jump points of initial data, the atom points of interaction 
and cancelation measure /i IC and the points where two sub-discontinuity curves of different families cross 
each other. 

Consider a point P = (r, £) ^ 0. Since u(-,t) has bounded variation, there exist the limits 

u L := lim u(x,t), u r := lim u{x, r). 

x— >{— x— >f+ 

Assuming that u K = Ti[u L ](s). We only consider the case for s > and the case for s < is analogous 
to prove. 

Applying the tame oscillation condition (see p. 295 of [3]), one obtains 
(5.1) lim u(x.t) — u L , lim u(x,t) — u R . 

(x,t)^( 5 , T ) (x,t)-H 5 ,x) 
x<t<T + (e-x)/A T<t<T + (x-0/X 

for some constant A which is larger than all characteristic speeds. 
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Suppose that there are i-subdiscontinuity curves y^ 1 , • • • , y\ l satisfying j/ . J € ^ J (e^) with ej > 0, 1 < 

./• 

2/i fel W = ---=2/ ; fci (r)=e 

and there is no other subdiscontinuity curve passing through the point P. This can be done because of 
the conclusion of Lemma 4.8 and the fact that P ^ 0. It is easy to show that ej must have the same 
sign, i.e. e jl e j2 > for any ji, j 2 € {1, • • • , Z}. 
By rearranging the index, we can assume that 

and 

(H) there is no y%° £ such that j/q° -< y^ 1 with (t ) = y*» (t ) or yf> -< y fe ° with y^(t ) = yf< (to). 

Step 1. By the definition, there exist y^ v , • • • , y\ l v e & itV such that — » t/ -' , Vj € {1, • • • , Z}. We 
claim that 



(5.2a) lim limsup 



(5.2b) lim limsup 

r-s-0+ i^-voo 



sup \u v (x,t) — u L \ I = 0, 



sup |ztj/(a;, i) — u fl 

\ (x,t)e-B(P,r) / 



= 0, 



where B(P,r) is a ball centred at P with radius r. 

Indeed, if (5.2a) is not true , by the first limit in (5.1) and u v —¥ u pointwise a.e., there exist two 
sequences of points P v . Q v converging to P and P„, Q y on the left of y\ such that the segment P V Q V is 
space-like and 

u(P„) -> u L 

and 

- «!/(Qi/)| > £o- 

It is not restrictive to assume that the direction P V Q V towards y\\- 

Let Kj(P v Q v ) be the total wave strength of wave-fronts of j-th family which cross the segment P V Q V . 
Then, one has Kj(P v Q v ) > cq for some j <E {1, • • • , d}. We consider the following three cases. 

Case 1. If j > i, we take the maximal forward generalized j-characteristic \ + through P v and minimal 
generalized ^'-characteristic \ ~ through Q w . 

If \ + arL d x~ interact each other at 0„ before hitting y^ v . We consider the region T y bounded by 
PvQu, X + an( i X- Since no fronts can leave T u through x + an( l X - By (3.5) and (3.9), we obtain that 
the amount of interaction and cancellation (IV) for u u within the closure of T v remains uniformly 
positive as v — > oo. 

If x + interacts y\ x v at A v and x~ interacts y\ x v at B v , we consider the region T v bounded by 
PvQv, X + , X~ and y\\- Then either there exists a constant < c < 1 such that /i£ c (r„) > Cq6o 
or there exists a constant < c' Q ' < 1 such that fronts with total strength lager than CqEq hitting A V B V . 
By (3.5) and (3.9) we can determine that n„ (FV) > e uniformly. 

For both above two cases, T v is contained in a ball B(P, rv) with f„ — > as v — ^ oo, which implies 
that p} {{P}) > 0. This is against the assumption P ^ 0. 

Case 2. If j < i, we consider the minimal backward generalized j-characteristic through the point 
P v and the maximal backward generalized j-characteristic through the point Q u . Then by the similar 
argument for the case j > i, we get fi IC ({P}) > against the assumptions. 

Case 3. If j = i and for any j' ^ i, 1 < j' < d, A,v (P v Qu) — > as v — > oo. We claim that the 
maximum of the strengths of all fronts which cross P V Q U tends to zero when v — > oo. 

If it is not true, since h.ji(P v Q v ) — > for f ^ i, there must be fronts of i-th family across P V Q V with 
uniformly large strength We assume that up to a subsequence, their (i, fco)-substrength are uniformly 
large for some ko G {!,••• , h}., that is \s*°\ > e for some e > 0. 
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Then by Definition 4.3 there must be, for some eo > 0, (eo, i, fco)-approximate subdiscontinuity curves 
y\° v which contains the wave fronts (we use the same notation as their strength for convenience) and 
since y ° v are uniformly Lipschitz continuous curves, up to a subsequence, there is a Lipschitz continuous 
curves y* , such that 

and Uq{t) = £. By Definition 4.7 we obtain j/q° -< y^ 1 , which contradicts the assumption (H). 
So we can always choose Q' v , P' v £ P V Q V such that 

U U {Q' U ) -> u L 

and 

K{Q' V ) - u v {Pl)\ > c e , 

where < cq < 1 and u l ,(Q' IJ ), u v (P'J) locate in the same Af for some k. 

Since for f ^ i, A J v(P^Q t/ ) is arbitrary small when v is large enough and the strength of fronts 
belonging to z-th family is small, one has 

max min \u v [H v ) - Ri[u v (P^)]{s*)\ < 1. 



Which means the values of each u v along the segment PIQ' V remain arbitrary close to the i-rarefaction 
curve through u L . Then by the analogous argument in the proof of Lemma 4.8, one gets the contradiction 
H IC ({P}) > 0. Therefore, we conclude that (5.2) is true. And (5.2b) is similar to prove. 

Step 2. Define 8? = \J i 3^. If P ^ n Graph( t $ 7 ) and if u is not continuous at P, then there exist 
e > and P„, Q v — > P such that P V Q V is space like and 

u v (P v ) -> u(P), \u v {Q v ) - u(P)\ > e for all v. 

Up to subsequence, we consider th following two cases. 

1) there exists j ^ f, such that min{Aj ;(P V Q V ), Aj/(P„Q„)} > eo- This situation can be ruled out 
by the argument in Case 1 of Step 1. 

2) For some j € {1, • • • , j} and all v, Aj(P u Q v ) > e and for all j' ^ j, Kj>{P v Q u ) — > as v — > oo. 
Then one can use the argument in Case 3 of Step 1 to obtain the contradiction. 

Therefore we get the continuity of u outside Graph(^) (J 9. 
Then, we have proved the first part of the theorem. 

Step 3. We now establish the Rankinc-Hugoniot condition (1.5) for curves in ST. Let P = (to,Xo) E 
Graphi y 5 r ) \ 0, and write 

u L = lim u(x,to), u R — lim u(x,to). 

x^xq— x^xo-\- 

We consider two cases. 

Case 1. There is only one curve y € 3\ passing through point P. From (5.2), we know that the 
discontinuity [it L ,u R ] must be simple. Suppose that 3 y u — > y as v — > oo. By (3.3) and (3.4), we 
obtain 

a i (^ ) <)[^-<] = /(^)-/«), 
\y v - ai(v£,u*)\ < 2e u , 

where 

= lim u(t u ,x), u R — lim u(t u ,x), t u — > to as v °°- 

x-*y u (t u )- x-H/„(t„)+ 

Then by (5.2), one has for every e > there exists D(e) such that W > v, one has 
|i/ I/ -o- i (u L ,u R )| < |l/ v -<7<(u^u?)| + |o-t(t» L ,« R )-o-i(«5,«*)| <e. 

From (4.9) and the fact that the z-waves of Riemann problem constructed by Theorem 2.1 are Liu 
admissible, one deduces that 

y(t )(u R -u L ) = f(u R )-f(u L ), 

and 

y(to) <^(5,[ U L ](r), U L ), Vre [0,s\. 
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Case 2. If the discontinuity [it L ,it R ] is a composition of [uo,ui], [1*1,1x2], • • • , [it;,u;+i] where uo 

= u R , Uj = Ti[u h ](sj) and u L G Af , u R G A^ 2 *. Let 

fci = minjfc is even, fe > fc^}, 

fc p = max{fc is even, fe < k^} 

and k p = ki + 2(p — 1). One has p > I + 1. 

According to (5.2), there exist P V ,Q V — > P such that 



and the segment P,^ is space-like. 

We claim that there exist p subdiscontinuity curves y±, • • • , y p G S\ passing through point P, where 
y 3 G ^ +2j , j = !,-■■ , p - 1. 

In fact, let S^(P V Q U ) denote the maximal (i, fc)-substrength of all fronts across P V Q V . It is sufficient 
to show that, up to a subsequence, there is a constant C > such that 

S^(P~M~ V )>C, VfcG {fci,-- - ,fc p }. 

If not, there exists fep G {fci, • • • , fc p } such that S^{P V Q V ) ->0asv^ 00. Since p IC (P) = 0, we have 



Aj(P„Q„) ->• 0, as !/ -> 00, Vj 7^ i. 

By Lemma 4.8, we conclude that S^{P V Q V ) — > 0, as 1/ — > 00, Vfc odd. Since the (i, fco)-substrength of all 
wave fronts are arbitrary small, up to a subsequence, one can always find points P' V ,Q' V on P V Q V such 
that u v {P'^), u v (Q' v ) G A^° and all fronts are either admissible discontinuity with left and right value 
inside A^° or rarefaction fronts. Therefore by the analogous argument in Case 3 of Step 1, there exist a 
constant c > independent on v such that in a small neighborhood T v of P, one has (r„) > c. This 
contradicts with the assumption p IC (P) = 0. This concludes our claim. 

Moreover, by Lemma 4.8 and the equalities (5.2), there are exactly p subdiscontinuity curves passing 
through P, otherwise P must be an interaction point with /j, IC (P) > 0. 

fe k n 

Suppose u n+ i — Ti[u n ]{s n ) for some s n > 0, n G {0, • • • , /} and Ti[u n }(-) intersects Z i ni , ■ ■ • , Z i q 

fe fe 

at Un,i, m " ,u n ,q- Then the subdiscontinuity curves with substrength s^" 1 ,--- ,s i " q must conincide in 
neighborhood of time r. 

If in a neigborhood of P, yj and yj+i are not identical, by the similar argument for proving (5.2), one 
can show that (1.8) is true. 

Suppose ^i^ v 3 yj. v — > yj, j G {1,-- ■ ,p}- As we discussed above, it is not restrictive to assume that 
p = I + 1 and there is a neighborhood U(to) of to, such that Vi € U(to), 

vtAt) < ■ ■ ■ < Vp,v 

Then, using the similar argument in Step 1, besides (5.2) and (5.2b), one can also show that 



lim limsup 

r->0+ v^rac 



sup \u v {x,t) - Um-il \= m = 2, ■ ■ ■ ,p, 



1 

\(x,t)eB(P,r) 

( 



lim limsup 



sup ti^a;, t) — m„| ] =0, n = I, ■ ■ ■ ,1. 

\(l,t)eB(P,f) 

For notational convenience, we write uo = « L , it p = u R . Therefore, by the same argument in Case 1, 
we obtain that for n — 1, ■ ■ • ,p 

yn{h){Un - U n -l) = f(.U n ) - f(u n -l) 

and 

y n (to) < <Ti(u n ,Si[u L ](T)), Vr G [sj-i,Sj]. 
Adding them together, one finally obtain for m = 1, ■ ■ • ,p 

ym{t ){u R -U L )=f{ U R )-f{u L ) 
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and 

y m (t ) < a l (u L ,S l [u L }{T)), Vr e [o,s]. 

□ 

The proof of Theorem 1.1 already implies the results of Theorem 1.2. 

Concerning a sequence of exact solution of (1.1) such that u v — > u in Lj , one can approximate each 
u v by a sequence of wave-front tracking approximations u mM —¥ u v and by a suitable diagonal sequence 
u u,m{u) ~~ ^ u one has the following corollary. 

Corollary 5.2 (stability of discontinuity curves). Considering a sequence of exact solutions u v such that 
u v — > u in L} , one has 

(1) Let y v : [t~,t+] t— > R be a discontinuity curves of u v described in Theorem 1.1. Assume t~ — > 
t~,t£ — > t + and y u {t) — > y(t) /or eac/i t G TTien y(-) is a discontinuity curve of the 
limiting solution u with the properties mentioned in Theorem 1.1. 

(2) Viceversa, let y : [t~,t+] i— > K fee a discontinuity curve of u for a.e. t £ TTien i/iere 
exists a sequence of discontinuity curves y v : [i~ , t^T] i— > K such that 

t~ -t t~, t+ -> t + , lim y v (t) = y(t), 

for a.e. t G [t~ , t + \. 

6. An remark on general strict hyperbolic systems 

We construct a strict hyperbolic system of conservation laws with one characteristic family which is 
not linearly degenerate or piecewise genuinely nonlinear. Therefore neither the assumption of Theorem 
1.1 or that of Theorem 10.4 in [6] holds. We show that the set of jump points of its admissible solution 
to some initial data can not be "exactly" covered by countably many Lipschitz continuous curves. 

Consider the following 2x2 system 



(6.1) 



\u t + f(u,v) x = 0, 
\v t -v x = 0. 

where / is a smooth function and u, v is the unknown variables. The Jacobian matrix of flux function is 

DF{u,v)=^ [\ 
Then the eigenvalues are 

Ai = — 1, A2 = /„. 
And the corresponding right eigenvectors are 

ri(«,«) = (/„,-/„ -If, r 2 = (l,0) T . 

The system is strict hyperbolic if f u > — 1. (In fact, / constructed latter satisfies this property). 
Obviously, one has 

Z x = {{u,v) : VAi • ri(u,v) = 0} = M 2 , 
which means that the 1-th characteristic family is linearly degenerate. 
Latter we will also show that 

(6.2) Z 2 = {(«,«) : VA 2 ■ r 2 } = {(«,«) : f uu (u,v) = 0} = {v - 0}. 

This yields that the vector field r 2 is tangent to the manifold Z2, therefore the second characteristic 
family is not piecewise genuinely nonlinear or linearly degenerate.. 

Define f(u,v) = e~ 1 ^ v u 2 /2 when v > and f(u,0) = 0. In the following, we define the value of / for 
v < 0. 

Let the initial data to be 



(6.3) u (x) = < v (x) 

\u r X > 

for some small constants ui > u r and a, h > 0. 
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Since the second equation in (6.1) is a linear transport equation, one has 
(6.4) v(x,t) 



—a x + t < h, 
a x + t > h. 



Then one can solve the system (6.1) by regarding it as a scalar conservation laws of u 

u t + f(u,v) x = 

with discontinuous coefficient when v is a piecewise constant function. 

If w.r.t. u, f is concave for some small fixed v < and concave for some fixed v > 0, then u is a center 
rarefaction wave in the area {x + t < h}. Immediately after the rarefaction wave crosses the characteristic 
line x + 1 = h, it turns out to be a compressive wave since f(u, —a) is a convex function on u. It may 
generate a shock after a short time. 

Indeed, this can be done in the following way. 

Consider a center rarefaction wave of u in the area {x + t < h} 

{ui x/t < f u (ui,-a), 

g(x/t) x/t = f u (g(x/t),-a), 
u r x/t > f u (u r ,-a), 

where <?(•) is the converse function of /«(•, —a). 

Let u~ = <?(£*) for some £*. The value of u will be u~ along the characteristic line 

x = tf u (u~,-a) 

until it intersects another characteristic line x + t = h. 
Solving the equations 



(6.6) 



\x = tf u (u ,-a), 
}x = —t + h, 

we get the intersection point of two characteristics: 

\xq = hf u (u~ ,-a)/[l + fu(u~ ,-a)] 



(6.7) 



t Q = h/[l + f u (u ,-a)}. 



Next, we compute the value of u after the characteristic cross the line x + t = h. For each point (xq, to) 
on the line x + t — h, let 

u + (xo,to) = lim u(x,t). 

x+t>h 
(x,t)-t(x ,t ) 

By Rankine-Hugoniot condition, one has 

-l/a( +\2 

-(«+ -u-) = - f(u-,-a). 

This yields 

(6.8) u+ = ^Wl + 2e-yHf(u-,-a) + u-)^ 

g— l/a 

bmce f(u,v) = e- 1 ' a u 2 /2 on area {h < x + 1 < 2h} , the characteristic line of the equation 

ut + f(u,a) x = 

starting at (xo,to) is 

x - h ff- L - a \ = e-V.„ +(t - " ). 

! + /«(« ,-a) 1 + f u \u ,-a) 

We require that it passes through the point (0, 2h), that is 

hf u (u-,-a) _ ^_ 1/a _ +/ot /i 



- e- 1/a w+(2/i — -), 

! + /„(«-, -a) V l + / u (u-,-a) ; ' 
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which yields 

(6.9) fu(u-,-a) = 



_ e -l/a u + 



2e -i/a u + +1 - 
Substitute (6.8) into (6.9), we get 

1 - g{u~,-a) 



f u {u , -a) = 



2g(u , -a) - 1 ' 

where g{u~ , —a) = \J\ + 2e~ 1 / a (f(u~ , —a) + 

Now, we consider the Cauchy problem of an ODE with parameter a 

^F(u,a) = 2G(Fu,a)-l ' 



(6.10) 



F(0,a)=0, 



where G(F, u, a) = y/l + 2e- 1 / a (F + u). 

By the theory of the ODE, since G is smooth when (u, a, F) lies in a small neighborhood of the origin, 
there is a unique smooth solution F defined on some interval [—6, 6] (b > 0), smoothly depending on the 
parameter a. 

Therefore we defined / for v < small as 

f(u,v) = F(u,v). 

By our construction, j{u, v) is concave about u for a negative fixed v. In fact, from (6.10), one finds 

f = — < 
Juu 2 g+l ' 



-l/v 



since g u — — ^ ^- This concludes that / is concave with respect to u. 

It is not difficult to verify that all partial derivatives of / is continuous on {v — 0}. As / is smooth on 
{v 7^ 0}, it concludes that / is smooth and independent of h. 




Figure 3. 

By the above construction, as shown in Figure 3. The rarefaction wave turns out to be compressive 
in the area {h < x + t < 2h} which eventually generates a shock starting at the point (0, 2h) and 
propagating along the x = 0. Notice that there is another shock starting from the point (h, 0). However, 
we can modify the initial data a little to get rid of this shock. In fact, recalling the formula (6.8) and 
letting 

-1 + \/T+ 2e- 1 / a (f(u r , -a) + u r ) 
e -l/a 
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We can replace uo in the initial data by 

{ui x < 0, 
u r < x < h, 
U\ x > h 

to make the speed of the characteristics starting from {h < x < 2h} the same as nearby. By the total 
variation estimates for the general system 

Tot.Var.{u(-,i)} < Tot.Var.{u (-)}, 

it is not restrict to assume that the total variation of uq is sufficiently small. 

Now, we begin to find the initial data with which (6.1) had an admissible solution containing " Cantor 
set" shocks. 

Now, we consider the initial data 

{0 x < 2h, 
-a 2h< x < 3h, 
<ih^ / Ah 
a on < x < An, 
x> Ah. 

Since the second equation in (6.1) is a linear transport equation, one has 

'0 x + t<2h, 

,(,1.2. r,r. !} = {- a ^<X + t<3h, 

' a 3h<x + t<Ah, 

x + 1 > Ah. 

In the area {x + 1 < 2h}, the first equation in (6.1) becomes ut = 0, then one has 

ui x < min{2/i — t, 0}, 



(6.13) u{x,t) 



u r x > max{2/i — t, 0} 



Next, we compute the value of u in {2h < x + t < 3h}. For any (x 1 , tf) ^ (0, 2h) on the line x + t = 2h, 
write 

ur = lim u(x,t), u + = lim u(x,t). 

(x,t)->(a:',i') (x,t)-¥(x' ,t') 

x+t<2h x+t>2h 

By Rankine-Hugoniot conditions, one has 

/(u+ -a)- /(u-Q) = -(«-- u~). 

which yields 

u~ = u + + f(u + , —a). 

Regarding u + as a function of u~ and differentiating the above equation on both sides with respect to 
u~, one gets 

(/„ + l)(u+) / = l. 

By (6.9), one has 

, 4-M 1 20- 1 

i + fu g 

in a small neighborhood of origin. Thus u + is strictly increasing w.r.t. u~ . This concludes that at point 
(0, 2h) the left value of u is still larger than the right value of u, i.e. 

where 



= lim u(x,t), = lim u(x,t). 



(x,t)->(0,2h) (x,t)^(0,2h) 
2h~t<x<0 i>max{2ft-i,0} 

As we discussed before, one gets a rarefaction wave by solving the Riemann problem in {2h < x + t < 3/i} 
and a compressive waves in area {3/i < x + 1 < Ah}. 
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Next, we compute the value of u in the zone {x + 1 > Ah}. For any point (xq, to) ^ (0, Ah) on the line 
x + 1 = Ah, let 

u + (xo,to)= Hm u(x,t), u~(xo,t )= lim u(x,t). 

x+t>4h x+t<4h 
(x.t)^(x ,t ) (x,t)^(x ,t ) 

By Rankine-Hugoniot conditions, one has 

f(u+,0)- f(u-,a) = -(u+ -u-). 
Then since f(u + ,0) = and f(u~,a) — e _1 / a (u~) 2 /2, one obtains 

, e-^^u-) 2 



We claim that u\ > u^, where 



u \ = lim u(x, t), 



O,t)->(4h,0) 
4h-t<x<0 



lim u(x,t). 

(x,t)^(4hfi) 
l>max{4h-t,0} 



are the left /right of u across the line x + t = Ah In fact, since u L > u Rl where 

u~[ = lim u(x, t) 



(x,t)^(4hfi) 

x+t<4h 
x<\(t-4h) 



u„ = lim u(x, t). 

(x,t)~f(4h,Q) 

x+t<4h 
x>-\(t-4h) 



one has 



U L - Uj 
,-1/a 



i- 1/a (u- R ) 2 



(u L + u R ) > 



e" 1/a («Z) 2 



since a and (u~[,u R ) sufficiently small. 

Therefore as the first equation in (6.1) is u t — in {a; + t > Ah}, the jump prorogate along x = 
which turns out to be a shock. Similarly, by modifying the initial data a little, we can guarantee that 
there is not other shocks in the solution. (See Figure 4.) 



rarefaction wave 



compressive wave 




2h 4h 



Figure 4. Afer modifying a little the initial data (6.11), the admissible solution contains 
two shocks along the i-axis. 
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Next, we construct an admissible solution to such that it consists of graph of shock curves which is an 
1-dimcnsional Cantor set. Let 

B-^i^^ 2 )-^, i£ n =(*^,^y 6l n = 0,l,-,3— -1. 
V 3™ 1 3 m / V 3 m 3 m I 

and 

3 m-i_ 1 

B m ,n — B mn U B^ l n and _B m = B m ^ n . 

n=0 

We consider initial datum as 

€ B mn , 

u , m (x) = { 1 2 ^ «o,m = "( -ct„ a: G £?+ n , 

a: G R\-B m ,„, 

where one can always choose ui, u r and a suitable sequence {a n } such that total variation of (uq.vq) 
is sufficiently small. By modifying uq properly similarly to get rid of extra shocks. One finds that the 
admissible solution (u m ,v m ) of the system (6.1) is continuous inside K rn where 

l°°l 3 T 1 |~ 1 r/ \ + 6(3n + l) 6(3n + 2) 

K m := |J U {(*>*) € x : 3m ~ * < * < 3TO - 4 

m— 1 n— 

and its shocks belongs to 1-th family are located on 

{(0,0: i£ |0,o O ,}\ I U>0,0:&^<<<&±l)}. 

fe=0 ) 

Obviously, (u,v) := lim m _> 00 (u m ,i' m ) is the admissible solution of (6.1) with initial data (uq,Vq) := 
lim(ito, m , fo,m)- Asm-) oo, C m := [0,6] \ B m converge to the Cantor set C (scaled by 6). Since the 
Cantor set is uncountable set and can not contain any interval of non-zero length, it is impossible to find 
countable many Lipschitz continuous curves which exactly cover discontinuities of (u, v), that is, on these 
curves, either (u, v) is discontinuous or it is an interaction point. This means that Theorem 1.1 fails for 
this situation. 
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